Towards an experimental von Karman dynamo: numerical studies for an optimized 

design 
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Numerical studies of a kinematic dynamo based on von Karman type flows between two counter- 
rotating disks in a finite cylinder are reported. The flow has been optimized using a water model 
experiment, varying the driving impellers configuration. A solution leading to dynamo action for 
the mean flow has been found. This solution may be achieved in VKS2, the new sodium experiment 
to be performed in Cadarache, France. The optimization process is described and discussed, then 
the effects of adding a stationary conducting layer around the flow on the threshold, on the shape 
of the neutral mode and on the magnetic energy balance are studied. Finally, the possible processes 
involved into kinematic dynamo action in a von Karman flow are reviewed and discussed. Among 
the possible processes we highlight the joint effect of the boundary-layer radial velocity shear and 
of the Ohmic dissipation localized at the flow/outer-shell boundary. 

PACS numbers: 47.65-|-a, 91.25.Cw 



I. INTRODUCTION 

In an electrically conducting fluid, kinetic energy can 
be converted into magnetic energy, if the flow is both of 
adequate topology and sufficient strength. This problem 
is known as the dynamo problem (jy, and is a magnetic 
seed-field instability. The equation describing the behav- 
ior of the magnetic induction field B in a fluid of resis- 
tivity ?7 under the action of a velocity field v writes in a 
dimensionless form: 



-^Vx(vxB) 



V*C* 



-V^B 



(1) 



where C* is a typical length scale and V* a typical 
velocity scale. In addition, one must take into account 
the divergence- free of B, the electromagnetic boundary 
conditions and the Navier-Stokes equations governing the 
fluid motion, including the back-reaction of the magnetic 
field on the flow through the Lorentz force. 

The magnetic Reynolds number i?„i — V*C*ri^^, 
which compares the advection to the Ohmic diffusion, 
controls the instability. Although this problem is sim- 
ple to set, it is still open. While some flows lead to the 
dynamo instability with a certain threshold i?^, other 
flows do not, and anti-dynamo theorems are not suffi- 
cient to explain this sensitivity to flow geometry 0] . The 
two recent experimental success of Karlsruhe and Riga 
3, 4, 5, 6J are in good agreement with analytical and 
numerical calculations 0j H 0, ^| ; these two dynamos 
belong to the category of constrained dynamos: the flow 
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is forced in pipes and the level of turbulence remains low. 
However, the saturation mechanisms of a dynamo are not 
well known, and the role of turbulence on this instability 
remains misunderstood [H El El El El El 113 • 

The next generation of experimental homogeneous un- 
constrained dynamos (still in progress, see for example 
Frick et al, Shew et al, Marie et al. and O'Connell et 
al. in the Cargese 2000 workshop proceedings (Tsll) might 
provide answers to these questions. The VKS liquid- 
sodium experiment held in Cadarache, France jl9ll2(ll21| 
belong to this category. The VKS experiment is based 
on a class of flows called von Karman type flows. In a 
closed cylinder, the fluid is inertially set into motion by 
two coaxial counterrotating impellers fltted with blades. 
This paper being devoted to the hydrodynamical and 
magnetohydrodynamical properties of the mean flow, let 
us first describe briefiy the phenomenology of such mean 
flow. Each impeller acts as a centrifugal pump: the 
fluid rotates with the impeller and is expelled radially 
by centrifugal effect. To ensure mass conservation the 
ffuid is pumped in the center of the impeller and recircu- 
lates near the cylinder wall. In the exact counterrotating 
regime, the mean flow is divided into two toric cells sep- 
arated by an azimuthal shear layer. Such a mean ffow 
has the following features, known to favor dynamo ac- 
tion: differential rotation, lack of mirror symmetry and 
presence of a hyperbolic stagnation point in the center of 
the volume. In the VKS experimental devices, the flow, 
inertially driven at kinetic Reynolds number up to 10^ 
(see below), is highly turbulent. As far as full numer- 
ical MHD treatment of realistic inertially driven high- 
Reynolds- number flows cannot be carried out, this study 
is restricted to the kinematic dynamo capability of von 
Karman mean flows. 

Several measurements of induced flelds have been per- 
formed in the first VKS device (VKSl) '2^, in rather 
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good agreement with previous numerical studies |22| , but 
no dynamo was seen: in fact the achievable magnetic 
Reynolds number in the VKSl experiment remained be- 
low the threshold calculated by Marie et al. A 
larger device — VKS2, diameter 0.6 m and 300 kW power 
supply — is under construction. The main generic prop- 
erties of mean flow dynamo action have been highlighted 
by Marie et al. p3 | on two different experimental von 
Karman velocity fields. Furthermore, various numeri- 
cal studies in comparable spherical flows confirmed the 
strong effect of flow topology on dynamo action . 
In the experimental approach, lots of parameters can be 
varied, such as the impellers blade design, in order to 
modify the flow features. In addition, following BuUard 
& Gubbins several studies suggest to add a layer 
of stationary conductor around the flow to help the dy- 
namo action. All these considerations lead us to consider 
the implementation of a static conducting layer in the 
VKS2 device and to perform a careful optimization of 
the mean velocity field by a kinematic approach of the 
dynamo problem. 

Looking further towards the real VKS2 experiment, 
one should discuss the major remaining physical unex- 
plored feature: the role of hydrodynamical turbulence. 
Such turbulence in an inertially-driven closed flow will 
be very far from homogeneity and isotropy. The pres- 
ence of hydrodynamical small scale turbulence could act 
in two different ways: on the one hand, it may increase 
the effective magnetic diffusivity, inhibiting the dynamo 
action l^^. On the other hand, it could help the dynamo 
through small-scale a-effect 26j. Moreover, the presence 
of a turbulent mixing layer between the two counterrotat- 
ing cells may move the instantaneous velocity field away 
from the time-averaged velocity field for large time-scales 
[27| . As the VKS2 experiment is designed to operate 
above the predicted kinematic threshold presented in this 
paper, it is expected to give an experimental answer to 
the question about the role of turbulence on the instabil- 
ity. Furthermore, if it exhibits dynamo, it will allow to 
study the dynamical saturation regime which is outside 
the present paper scope. 

In this article, we report the optimization of the time- 
averaged fiow in a von Karman liquid sodium experi- 
ment. We design a solution which can be experimentally 
achieved in VKS2, the new device held in Cadarache, 
France. This solution particularly relies on the addi- 
tion of a static conducting layer surrounding the flow. 
The paper is organized as follows. In Section [H] we first 
present the experimental and numerical techniques that 
have been used. In Section fllll we present an overview of 
the optimization process which lead to the experimental 
configuration chosen for the VKS2 device. We study the 
infiuence of the shape of the impellers both on the hy- 
drodynamical fiow properties and on the onset of kine- 
matic dynamo action. In Section IIVI we focus on the 
understanding of the observed kinematic dynamo on a 
magnetohydrodynamical point of view: we examine the 
structure of the eigenmode and the effects of an outer 



conducting boundary. Finally, in Section El "we review 
some possible mechanisms leading to kinematic dynamo 
action in a von Karman flow and propose some conjec- 
tural explanations based on our observations. 

II. EXPERIMENTAL AND NUMERICAL 
TOOLS 

A. What can be done numerically 

The bearing of numerical simulations in the design of 
experimental fluid dynamos deserves some general com- 
ments. Kinetic Reynolds numbers of such liquid sodium 
flows are typically 10^, well beyond any conceivable di- 
rect numerical simulation. Moreover, to describe effec- 
tive MHD features, it would be necessary to manage very 
small magnetic Prandtl numbers, close to 10~^, a value 
presently out of computational feasibility. Several groups 
are progressing in this way on model flows, for exam- 
ple with Large Eddy Simulations which can reach 
magnetic Prandtl numbers as low as 10~^- 10~^. An- 
other strong difficulty arises from the search of realistic 
magnetic boundary conditions treatment which prove in 
practice also to be difficult to implement, except for the 
spherical geometry. 

An alternative numerical approach is to introduce a 
given flow in the magnetic induction equation ^ and to 
perform "kinematic dynamo" computations. This fiow 
can be either analytical H, |2^ , computed by pure hy- 
drodynamical simulations (which may now be performed 
with Reynolds numbers up to a few thousands), or 
measured in laboratory water models (22k .28,] by Laser 
Doppler velocimetry (LDV) or by Particle Imaging Ve- 
locimetry (PIV). Such measurements lead to a map of 
the time-averaged fiow and to main properties of the fiuc- 
tuating components: turbulence level, correlation times, 
etc... Kinematic dynamo computations have been suc- 
cessfully used to describe or optimize the Riga and 
Karlsruhe 8] dynamo experiments. 

We will follow here the kinematic approach using the 
time-averaged flow measured in a water model at real- 
istic kinetic Reynolds number. Indeed, potentially im- 
portant features such as velocity fluctuations will not be 
considered. Another strong limitation of our pragmatic 
kinematic approach is its linearity: computations may 
predict if an initial seed field grows, but the study of the 
saturation regime will rely exclusively on the results of 
the real MHD experiment VKS. 

B. Experimental measurements 

In order to measure the time-averaged velocity field — 
hereafter simply denoted mean field — we use a water- 
model experiment which is a half-scale model of the 
VKS2 sodium device. The experimental setup, measure- 
ment techniques, and methods are presented in detail in 



3 



Na at rest 



Liquid Na 



IS 



w 



Hr 



FIG. 1: Sketch of the VKS2 experiment. The container radius 
Rc is taken as unit scale, w is the dimensionless thickness of 
sodium at rest. 



Refs. 123,123. However, we present below an overview 
of our experimental issues and highlight the evolutions 
with respect to those previous works. 

We use water as working fluid for our study, noting 
that its hydrodynamical properties at 50°C (kinematic 
viscosity v and density p) are very close to sodium prop- 
erties at 120°C. 

A sketch of the von Karman experiments is presented 
in Fig. n The cylinder is of radius Rc and height 
He — l.SRc- In the following, all the spatial quanti- 
ties are given in units of Rc = C* . The hydrodynami- 
cal time scale is based on the impeller driving frequency 
/: if V is the measured velocity field for a driving fre- 
quency /, the dimensionless mean velocity field is thus 
V = (2^i?,/)-iV. 

The integral kinetic Reynolds number Re is typically 
10^ in the water-model, and 10'' in the sodium device 
VKS2. The inertially driven flow is highly turbulent, 
with velocity fluctuations up to 40 percent of the maxi- 
mum velocity |2(i .22J . In the water model, we measure 
the time-averaged velocity field by Laser Doppler Ve- 
locimetry (LDV). Data are averaged over typically 300 
disk rotation periods. We have performed measurements 
of velocity in several points for several driving frequen- 
cies: as expected for so highly turbulent a flow, the di- 
mensionless velocity v does not depend on the integral 
Reynolds number Re = V*C*v-^ [13 • 

Velocity modulations at the blade frequency have been 
observed only in and very close to the inter-blade do- 
mains. These modulations are thus time-averaged and 
we can consider the mean flow as a solenoidal axisym- 
metric vector field |33| . So the toroidal part of the veloc- 
ity field Ve (in cylindrical coordinates) and the poloidal 
part {Vz^Vr) are independent. 

In the water model experiment dedicated to the study 
reported in this paper, special care has been given to 
the measurements of velocity fields, especially near the 
blades and at the cylinder wall, where the measurement 
grid has been refined. The mechanical quality of the 
experimental setup ensures good symmetry of the mean 
velocity fields with respect to rotation of tt around any 
diameter passing through the center of the cylinder {TZt,- 




FIG. 2: Dimensionless mean velocity field measured by LDV 
and symmetrized for kinematic dynamo simulations. Cylin- 
der axis is horizontal. Arrows correspond to poloidal part 
of the flow, color code to toroidal part. We use cylindrical 
coordinates (r,9,z), with origin at the center of the cylinder. 



symmetry). The fields presented in this paper are thus 
symmetrized by TZ-,^ with no noticeable changes in the 
profiles but with a slightly improved spatial signal-to- 
noise ratio. With respect to Ref. |[22|, the velocity fields 
are neither smoothed, nor stretched to different aspect 
ratios. 

Fig. 12 shows the mean flow produced by the optimal 
impeller. The mean flow respects the phenomenology 
given in the Introduction: it is made of two toroidal cells 
separated by a shear layer, and two poloidal recirculation 
cells. High velocities are measured in the whole volume: 
the inertial stirring is actually very efficient. Typically, 
the average over the flow volume of the mean velocity 
held is of order of 0.3 x (27ri?c/). 

In addition to velocity measurements, we perform 
global power consumption measurements: torques are 
measured through the current consumption in the mo- 
tors given by the servo drives and have been calibrated 
by calorimetry. 



C. Kinematic dynamo simulations 

Once we know the time-averaged velocity field, we in- 
tegrate the induction equation using an axially periodic 
kinematic dynamo code, written by J. Leorat |32l |. The 
code is pseudo-spectral in the axial and azimuthal direc- 
tions whether radial dependence is treated by high-order 
finite difference scheme. The numerical resolution cor- 
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responds to a grid of 48 points in the axial direction, 4 
points in the azimuthal direction (corresponding to wave 
numbers ni — 0, ±1) and 51 points in the radial direction 
for the flow domain. This spatial grid is the common ba- 
sis of our simulations and has been refined in some cases. 
The time scheme is second-order Adams-Bashforth with 
diffusive time unit td — R^r]~^. Typical time step is 
5.10~^ and simulations are generally carried out over 1 
time unit. 

Electrical conductivity and magnetic permeability are 
homogeneous and the external medium is insulating. Im- 
plementation of the magnetic boundary conditions for a 
finite cylinder is difficult, due to the non-local character 
of the continuity conditions at the boundary of the con- 
ducting fluid. On the contrary, axially periodic boundary 
conditions write down easily, since the harmonic external 
field has then an analytical expression. We choose thus 
to look for axially periodic solutions, using a relatively 
fast code, which allows to perform parametric studies. To 
validate our choice, we compared our results with results 
from a finite cylinder code (F. Stefani, private commu- 
nication) for some model flows and a few experimental 
flows. It happens that in all these cases, the periodic and 
the finite cylinder computations give comparable results. 
This remarkable agreement may be due to the peculiar 
flow and to the magnetic eigenmodes symmetries: we 
do not claim that it may be generalized to other flow 
geometries. Indeed, the numerical elementary box con- 
sists of two mirror-symmetric experimental velocity fields 
in order to avoid strong velocity discontinuities along z 
axis. The magnetic eigenmode could be either symmetric 
or antisymmetric towards this artificial mirror symmetry 
181^ . In the quasi totality of our simulations, the mag- 
netic field is mirror-antisymmetric, and we verify that no 
axial currents cross the mirror boundary. The few ex- 
otic symmetric cases we encountered cannot be used for 
experiments optimization. 

Further details on the code can be found in Ref. [s^ . 
We use mirror-antisymmetric initial magnetic seed field 
optimized for fast transient p3 |. Finally, we can act 
on the electromagnetic boundary conditions by adding 
a layer of stationary conductor of dimensionless thick- 
ness w, surrounding the fiow exactly as in the experiment 
(Fig. This extension is made keeping the grid radial 
resolution constant (51 points in the flow region). 



measure the mean velocity field and the power consump- 
tion. Then we simulate the kinematic dynamo problem. 
We try to identify features favoring dynamo action and 
modify parameters in order to reduce the threshold and 
the power consumption and go back to the loop. 

B. Impeller tunable parameters. 

The impellers are fiat disks of radius R fitted with 8 
blades of height h. The blades are arcs of circle with a 
curvature radius C and are radial at the center of the 
disks. We use the angle a = arcsin(^) to label the dif- 
ferent curvatures (see Fig. O. For straight blades a = 0. 
By convention, we use positive values to label the di- 
rection corresponding to the case where the fiuid is set 
into motion by the convex face of the blades. In order 
to study the opposite curvature (a < 0) we just rotate 
the impeller in the other direction. The two counter- 
rotating impellers are separated by He, the height of 
the cylinder. We fixed the aspect ratio Hc/Rc of the 
fiow volume to 1.8 as in VKS device. In practice we 
successively examine the effects of each parameter h, R 
and a on some global quantities characterizing the mean 
fiow. We then varied the parameters one by one, until we 
found a relative optimum for the dynamo threshold. We 
tested 12 different impellers, named TMxx, with three 
radn (i? = 0.5,0.75 & 0.925), various curvature angles a 
and different blade heights h. 




R 



FIG. 3: Sketch of the impeller parameters. 7? is the dimen- 
sionless radius, a the blade curvature angle. The sign of a 
is determined by the sense of rotation: positive when rotated 
anticlockwise. 



III. OPTIMIZATION OF THE VKS 
EXPERIMENT 

A. Optimization process 

The goal of our optimization process is to find the im- 
peller whose mean velocity field leads to the lowest i?^ 
for the lowest power cost. We have to find a solution fea- 
sible in VKS2, i.e. with liquid sodium in a 0.6 m diameter 
cylinder with 300 kW power supply. We performed an 
iterative optimization loop: for a given configuration, we 



C. Global quantities and scaling relations 

We know from empirical results [2^ l2^ l28l | that the 
poloidal to toroidal ratio F of the flow has a great impact 
on the dynamo threshold. Moreover, a purely toroidal 
flow is unable of sustained dynamo action l34l l35l , while 
it is possible for a purely poloidal fiow |36l l37j. We also 
notice that, for a Ponomarenko flow, the pitch parameter 
plays a major role Jud^dli- All these results lead us to 
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first focus on the ratio 2.5 - 

r- 

where (P) is the spatially averaged value of the poloidal g _ 

part of the mean flow, and (T) the average of the toroidal q 
part. I 

Another quantity of interest is the velocity factor V: "I 
the dimensionless maximum value of the velocity. In our 
simulations, the magnetic Reynolds number is based g g _ 

on the velocity factor, i.e. on a typical measured velocity 
in order to take into account the stirring efficiency: 



^ max(||V||: 

2 TT f 



-^0 



T 
T 



-45 45 90 

a 



Rm = 27rR2fV/77 

We also define a power coefficient Kp by dimensional 
analysis. We write the power V given by a motor to 
sustain the flow as follows: 

P = Kp{Re, geometry)p R^ 

with p the density of the fluid and Q = 27r/ the driving 
pulsation. We have checked 29] that Kp does not depend 
on the Reynolds number Re as expected for so highly 
turbulent inertially driven flows (30l |. 

The velocity factor measures the stirring efficiency: the 
greater V, the lesser the rotation frequency needed to 
reach a given velocity. Besides, the lesser Kp, the lesser 
the power to sustain a given driving frequency. The di- 
mensionless number which we need to focus on compares 
the velocity effectively reached in the flow to the power 
consumption. We call it the MaDo number: 

V 

MaDo = — — 

The greater MaDo, the lesser the power needed to reach 
a given velocity (i.e. a given magnetic Reynolds num- 
ber). The MaDo number is thus a hydrodynamical effi- 
ciency coefficient. To make the VKS experiment feasible 
at laboratory scale, it is necessary both to have great 
MaDo numbers and low critical magnetic Reynolds num- 
bers i?^. The question laying under the process of opti- 
mization is to know if we could on the one hand find a 
class of impellers with mean flows exhibiting dynamo ac- 
tion, and on the other hand if we could increase the ratio 
MaDo/i?^. It means that we have to look both at the 
global hydrodynamical quantities and at the magnetic 
induction stability when varying the impellers tunable 
parameters h, R and a. 

Fig. 21 presents MaDo for the whole set of impellers. 
For our class of impellers, the MaDo number remains of 
the same order of magnitude within ±10%. Only the 
smallest diameter impeller (R = 0.5) exhibits a slightly 
higher value. In the ideal case of homogeneous isotropic 



FIG. 4: MaDo number vs a for all the impellers we have 
tested. R = 0.925(T), R = 0.75(B) and R = 0.5(«). Closed 
symbols: h — 0.2. Open symbols: h < 0.1 

turbulence, far from boundaries, we can show that what 
we call the MaDo number is related to the Kolmogorov 
constant Ck — 1.5 j38||. The Kolmogorov constant is 
related to the kinetic energy spatial spectrum: 

E{k) = Ck e'/' fc"'/' 

where e is the massic dissipated power, and k the wave 
number. If we assume that e is homogeneous — V being 
the total dissipated power we measure — we have: 

- ^ 
PttR^.H, 

Using the definition 




E{k)dk 



and assuming ^(f^) — and using the steepness of 
the spectrum, we obtain: 

Eik^) = iv2fc-i 

with fco — 27r/i?c the injection scale. Then the relation 
between the MaDo number and Ck would be: 

MaDo^ ~ 3^-4/3(^^^)-2/3Ck ^ 0.44Ck 
Rc 

i.e., with Ck = 1-5, we should have, for homogeneous 
isotropic turbulence MaDo ~ 0.81. In our closed system 
with blades, we recover the same order of magnitude, 
and the fact that MaDo does merely not depend on the 
driving system. Thus, there is no obvious optimum for 
the hydrodynamical efficiency. Between various impellers 
producing dynamo action, the choice will be dominated 
by the value of the threshold i?^^. 
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Let us first get rid of the effect of the blade height 
h. The power factor Kp varies quasi-linearly with h. 
As MaDo is almost constant, smaller h impellers require 
higher rotation frequencies, rising technical difficulties. 
We choose h = 0.2, a compromise between stirring ef- 
ficiency and the necessity to keep the free volume suffi- 
ciently large. 

D. Influence of the poloidal/toroidal ratio F 

In our cylindrical von Karman flow without conducting 
layer {w = 0), there seems to be an optimal value for F 
close to 0.7. Since the mean flow is axisymmetric and 
divergence-free, the ratio F can be changed numerically 
by introducing an arbitrary multiplicative factor on, say, 
the toroidal part of the velocity field. In the following, 
Fq stands for the experimental ratio for the measured 
mean velocity field, whereas F stands for the numerically 
adjusted ratio. 
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FIG. 5: Magnetic energy growth rate a vs. numerical ratio F. 
Rm ~ 100, w — 0. Simulations performed for two different 
mean velocity fields (impellers TM71 (A) and TM73 (T) of 
radius R = 0.75). Larger symbols correspond to natural Fq of 
the impeller. Vertical dashed line corresponds to optimal F — 
0.7. Closed symbols stand for stationary regimes, whereas 
open symbols stand for oscillating regimes for F < 0.6. 

In Fig. [3 we plot the magnetic energy growth rate a 
(twice the magnetic field growth rate) for different val- 
ues of F, for magnetic Reynolds number R„i = 100 and 
without conducting layer (w — 0). The two curves cor- 
respond to two different mean velocity fields which have 
been experimentally measured in the water model (they 
correspond to the TM71 and TM73 impellers, see table 
n for their characteristics). We notice that the curves 
show the same bell shape with maximum growth rate at 
F ~ 0.7, which confirms the results of Ref. p^ . 

For F < 0.6, oscillating damped regimes (open symbols 
in Fig. [SJ are observed. We plot the temporal evolution 



of the magnetic energy in corresponding case in Fig. 
these regimes are qualitatively different from the oscillat- 
ing regimes already found in 22] for non 7?.7r-symmetric 
F — 0.7 velocity fields, consisting of one mode with a 
complex growth rate: the magnetic field is a single trav- 
eling wave, and the magnetic energy, integrated over the 
volume, evolves monotonically in time. 
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FIG. 6: Typical damped oscillating regime for impeller TM70 
at F = 0.5, w = 0, Rm = 140. (a): temporal evolution of the 
magnetic energy E = J B'^. Straight line is a linear fit of the 
form E{t) = Eo expiat) and gives the temporal growth rate 
(7 = —12.1. (b): temporal evolution of the z component of 
B at the point r = 0.4, 6 = 0,z = —0.23 with a nonlinear 
fit of the form: Bz{t) — a exp(at/2) cos{Lut + <j)) which gives 
a = -12.2 and lo = 20.7. 

In our case, the velocity field is 7?,7r-symmetric. It is 
known that symmetries strongly constrain the nature of 
eigenvalues and eigenmodes of linear stability problems. 
Due to 7?,T-symmetry invariance of the evolution operator 
for the magnetic field, two type of solutions are allowed 

M 

• One 7?,7r-symmetric eigenmode with a real eigen- 
value. The corresponding bifurcation is steady. 
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• Two eigenmodes images one of the other by 7?.7r, 
associated with complex-conjugate eigenvalues. 

For r > 0.6, we always observed stationary regimes. Oth- 
erwise, for r < 0.6, starting the temporal integration 
with an initial condition for the magnetic field which has 
non vanishing projection on both eigenmodes, we obtain 
a mix of two modes with complex-conjugate growth rates 
and the magnetic energy decays exponentially while pul- 
sating (Fig. El). The same feature has been reported for 
analytical "s2^2 ~ flows" in a cylindrical geometry 
with a Galerkin analysis of neutral modes and eigenvalues 
for the induction equation 3£j. A major interest of the 
latter method is that it gives the structure of the modes: 
one mode is localized near one impeller and rotates with 
it, the other being localized and rotating with the other 
impeller. Growing oscillating dynamos are rare in our 
system: a single case has been observed, for TM71(— ) 
(Fo = 0.53) with aw = 0.4 conducting layer at Rm ~ 215 
— 197, see table OJ. Such high a value for the mag- 
netic Reynolds number is out of the scope of our experi- 
mental study, and is close to the practical upper limit of 
the numerical code. 

Experimental dynamo action will thus be searched in 
the stationary regimes domain F > 0.6. Without con- 
ducting layer, we have to look for the optimal impeller 
around Fq ~ 0.7. 



E. Effects of the impeller radius R 
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One could a priori expect that a very large impeller 
is favorable to the hydrodynamical efficiency. This is 
not the case. For impellers with straight blades, MaDo 
slightly decreases with R: for respectively R = 0.5, 0.75 
and 0.925, we respectively get MaDo = 2.13, 1.64 and 
1.62. This tendency is below the experimental error. We 
thus consider that MaDo does not depend on the im- 
peller. 

Nevertheless one should not forget that V varies quasi- 
linearly with impeller radius R: if the impeller becomes 
smaller it must rotate faster to achieve a given value for 
the magnetic Reynolds number, which may again cause 
mechanical difficulties. We do not explore radii R smaller 
than 0.5. 

Concerning the topology of the mean flow, there is no 
noticeable effects of the radius R on the poloidal part. 
We always have two toric cells of recirculation, centered 
at a radius close to 0.75 ± 0.02 and almost constant 
for all impellers (see right part of Fig. [7||. The fluid is 
pumped to the impellers for < r < rp and is reinjected 
in the volume fp < r < 1. This can be interpreted as a 
geometrical constraint to ensure mass conservation: the 
circle of radius r — ^ (very close to 0.75) separates the 
unit disk into two regions of same area. 

The topology of the toroidal part of the mean flow now 
depends on the radius of the impeller. The radial profile 
of vg shows stronger departure from solid-body rotation 



FIG. 7: Radial profiles of toroidal velocity vg ((a)-(d)) for z — 
0.3 (dotted line), 0.675 (dashed line), & 0.9 (solid line); and 
axial velocity Vz ((e)-(h)) for various equidistant z between 
the two rotating disks. From top to bottom: experimental 
flow for (a-e): R = 0.5, (b-f): R = 0.75, (c-g): R = 0.925 
impeller and (d-h): model analytical flow given by equations 
(1) (see discussion below p. I15II . 

for smaller R (left part of Fig.[7|): this will be emphasized 
in the discussion. We performed simulations for three 
straight blades impellers of radii R = 0.5, R = 0.75 and 
R = 0.925; without conducting shell {w = 0) and with 
a conducting layer of thickness w — 0.4. We have inte- 
grated the induction equation for the three velocity fields 
numerically set to various F and compared the growth 
rates. The impeller of radius R — 0.75 close to the radius 
of the center of the poloidal recirculation cells systemat- 
ically gets greatest growth rate. So, this radius R = 0.75 
has been chosen for further investigations. 

F. Seek for the optimal blade curvature 

The hydrodynamical characteristics of the impellers of 
radius R = 0.75 are given in tabled For increasing blade 
curvature the average value of the poloidal velocity (P) 
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increases while the average value of the toroidal veloc- 
ity (T) decreases: the ratio Tq is a continuous growing 
function of curvature a (Fig. IHl. A phenomenological 
explanation for (T) variation can be given. The fluid 
pumped by the impeller is centrifugally expelled and is 
constrained to follow the blades. So, it exits the impeller 
with a velocity almost tangent to the blade exit angle 
a. Thus, for a < (resp. a > 0), the azimuthal ve- 
locity is bigger (resp. smaller) than the solid body rota- 
tion. Finally, it is possible to adjust Fg to a desired value 
by choosing the good curvature a, in order to lower the 
threshold for dynamo action. 




FIG. 8: Fo vs a for four impellers of radius R — 0.75 rotated 
in positive and negative direction (see Tabled. 

Without conducting shell, the optimal impeller is the 
TM71 (Fq = 0.69). But its threshold R^^ = 179 cannot 
be achieved in the VKS2 experiment. So, we now have to 
find another way to reduce i?^, the only relevant factor 
for the optimization. 
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FIG. 9: Shift in the optimal value of F when adding a con- 
ducting layer. Magnetic energy growth rate a vs. F for w — 
(•) and w = 0.4 (t). Impeller TM73, i?™ = 100. Larger 
symbols mark the natural Fo of the impeller. 
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G. Optimal configuration to be tested in the VKS2 
sodium experiment 

As in the Riga experiment |j, Ij, and as in numerical 
studies of various flows Eji > '^^ consider a sta- 
tionary layer of fluid sodium surrounding the flow. This 
significantly reduces the critical magnetic Reynolds num- 
ber, but also slightly shifts the optimal value for F. We 
have varied w between u; = and w — 1] since the ex- 
perimental VKS2 device is of fixed overall size (diameter 
0.6 m), the flow volume decreases while increasing the 
static layer thickness w. A compromise between this con- 
straint and the effects of increasing w has been found to 
he w — 0.4 and we mainly present here results concerning 
this value of w. In Fig. |21 we compare the bell-shaped 
curves obtained by numerical variation of the ratio F for 
the same impeller at the same in the case w — 0, 
and w = 0.4. The growth rates are much higher for 



FIG. 10: Growth rate a of magnetic energy vs numerical ratio 
F. Rrr, = 43, w = 0.4 for 4 different R = 0.75 impellers: TM70 
(•), TM71 (A), TM73 (T) and TM74 (►). Larger symbols 
mark the natural Fo of each impeller. 



w = 0.4, and the peak of the curve shifts from 0.7 to 0.8. 
We have performed simulations for four different velocity 
fields (Fig. Unjl, for w = 0.4 at i?„, = 43: the rcsuh is 
very robust, the four curves being very close. 

In Fig. ^] we plot the growth rates a of the magnetic 
energy simulated for four real mean velocity fields at var- 
ious Rm and for w — 0.4. The impeller TM73 was de- 
signed to create a mean velocity field with Fq = 0.80. It 
appears to be the best impeller, with a critical magnetic 
Reynolds number of i?^ = 43. Its threshold is divided 
by a factor 4 when adding a layer of stationary conduc- 
tor. This configuration (TM73, w = 0.4) wiU be the 
first one tested in the VKS2 experiment. The VKS2 ex- 
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TABLE L Global hydrodynamical dimensionless quantities (see text for definitions) for the radius 7? — 0.75 impeller family, 
rotating anticlockwise (+), or clockwise ( — ) (see Fig. |3J. The last two columns present the thresholds for kinematic dynamo 
action with {w = 0.4) and without (w — 0) conducting layer. Optimal values appear in bold font. Most negative curvatures 
have not been investigated (n.i.) but TM71( — ), which presents oscillatory (o) dynamo instability for R"^ = 197 with w = 0.4. 
(1): TM70 impeller (Fq = 0.60) has a tricky behavior exchanging stability between steady modes, oscillatory modes and a 
singular mode mirror-symmetric with respect to the periodization introduced along z and thus not physically relevant. 
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FIG. 11: Growth rate a vs natural ratio Fq for five impellers 
at various Rm and w = 0.4. From left to right: TM71(— ) 
with Fo = 0.53, TM70 (Fo = 0.60), TM71 (Fo = 0.69), TM73 
(Fo = 0.80), TM74 (Fo = 0.89), see also table|IJ. Closed sym- 
bols: stationary modes. Open symbols: oscillating modes. 



periment will be able to reach the threshold of kinematic 
dynamo action for the mean part of the flow. Meanwhile, 
turbulence level will be high and could lead to shift or 
even disappearance of the kinematic dynamo threshold. 
In Section IIVI we examine in details the effects of the 
boundary conditions on TM73 kinematic dynamo. 



address the question of the links between F and mean 
hehcity (H). In our case, for aspect ratio Hc/Rc — 1.8 
and impellers of radius R — 0.75, the mean helicity at 
a given rotation rate (H) = J v.(V x v) rdrdz does not 
depend on the blade curvature (see TableHJ. Observation 
of Fig. 1121 also reveals that the dominant contribution in 
the helicity scalar product is the product of the toroidal 
velocity {vg oc (T)) by the poloidal recirculation cells 
vorticity ((V x v)e cx {P))- We can therefore assume 
the scaling {H) oc {P){T), which is consistent with the 
fact that the product (P)(T) and (H) are both almost 
constant (Table QJ. 
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FIG. 12: Isocontour of kinetic helicity H - 
TM73 velocity field, (a): total helicity. (b): 
tribution ^^.(V x v)e is dominant. 



= v.(V X v) for 
azimuthal con- 



H. Role of flow helicity vs. Poloidal/Toroidal ratio 

Most l arg e scale dynamos known are based on helical 
flows As a concrete example, while successfully 

optimizing the Riga dynamo experiment, Stefani et al. [3 
noticed that the best flows were helicity maximizing. The 
first point we focused on during our optimization process, 
i.e., the existence of an optimal value for F, leads us to 



To compare the helicity content of different flows, we 
now consider the mean helicity at a given Rm, {H)/V^, 
more relevant for the dynamo problem. Figure 1131 
presents (i?)/V^ versus Tq for all h — 0.2 impellers. The 
R — 0.75 family reaches a maximum of order of 1 for 
Fq ~ 0.9. This tendency is confirmed by the solid curve 
which stands for a numerical variation of F for TM73 ve- 
locity field and is maximum for F = 1. Besides, even if 
R = 0.925 impellers give reasonably high values of helic- 
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ity near T = 0.5, there is an abrupt break in the tendency 
for high curvature: TM60 (see Ref. |22|) exhibits large 
To = 0.9 but less hehcity than TM74. Inset in Fig. HSl 
highlights this optimum for (H) / versus impeller ra- 
dius R. This confirms the impeller radius R — 0.75 we 
have chosen during the optimization described above. 
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FIG. 13: Mean helicity at a given Rn^ ((//>/V^) vs. poloidal 
over toroidal ratio. The R — 0.75 impeller serie (t) is plotted 
vs. To. The large open symbol stands for TM73 at To and the 
solid line stands for the same quantity plotted vs. numerical 
variation of TM73 velocity field (F). We also plot (H) /V^ vs. 
To for the R = 0.5 (★) and R = 0.925 (□) impellers. The 
inset presents (H) /V^ vs. impeller radius R for impellers of 
0.8 < To < 0.9. 



As far as the optimal value toward dynamo action for 
the ratio T (close to 0.7 — 0.8, depending on w) is lower 
than 1, the best velocity field is not absolutely helicity- 
maximizing. In other words, the best dynamo flow con- 
tains more toroidal velocity than the best helical flow. 
As shown by Leprovost |4l|, one can interpret the opti- 
mal r as a quantity that maximizes the product of mean 
helicity by a measure of the w-effect, i.e., the product 
{H){T) ^ {P){T)\ 



IV. IMPACT OF A CONDUCTING LAYER ON 
THE NEUTRAL MODE AND THE ENERGY 
BALANCE FOR THE VKS2 OPTIMIZED 
VELOCITY FIELD 

In this section, we deal with the mean velocity field 
produced between two counterrotating TM73 impellers 
in a cylinder of aspect ratio = 1.8, like the first exper- 
imental configuration chosen for the VKS2 experiment. 
See Table for the characteristics of this impeller, and 
Fig. 121 for a plot of the mean velocity field. We detail the 
effects of adding a static layer of conductor surrounding 
the flow and compare the neutral mode structures, the 
magnetic energy and current density spatial repartition 
for this kinematic dynamo. 



A. Neutral mode for it; = 

Without conducting layer, this flow exhibits dynamo 
action with a critical magnetic Reynolds number 7?^ = 
180. The neutral mode is stationary in time and has a 
TO = 1 azimuthal dependency. In Fig. 1141 we plot an 
isodensity surface of the magnetic energy (50% of the 
maximum) in the case w — at Rm — Rm — 180. The 
field concentrates near the axis into two twisted banana- 
shaped regions of strong axial field. Near the interface 
between the flow and the outer insulating medium, there 
are two small sheets located on both sides of the plane 
z = where the magnetic field is almost transverse to 
the external boundary and dipolar. The topology of the 
neutral mode is very close to those obtained by Marie et 
al. Il^l with different impellers, and to those obtained 
on analytical 52^2 ~ flows in a cylindrical geometry 
with the previously described Galerkin analysis [39j |. 

In Fig. El we present sections of the B and j fields, 
j = V X B being the dimensionless current density. The 
scale for B is chosen such as the magnetic energy in- 
tegrated over the volume is unity. Since the azimuthal 
dependency is m = 1, two cut planes are sufficient to 
describe the neutral mode. In the bulk where twisted- 
banana-shaped structures are identified, we note that 
the toroidal and poloidal parts of B are of the same or- 
der of magnitude and that B concentrates near the axis, 
where it experiences strong stretching due to the stagna- 
tion point in the velocity field. Around the center of the 




FIG. 14: Isodensity surface of magnetic energy (50% of the 
maximum) for the neutral mode without conducting layer 
(w = 0). Cylinder axis is horizontal. Arrows stand for the 
external dipolar field source regions. 
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FIG. 15: Meridian sections of B and j fields for the neutral mode with w = 0. B is divided by the total magnetic energy. 
Arrows correspond to components lying in the cut plane, and color code to the component transverse to the cut plane. A unit 
arrow is set into each figure lower left corner, (a): B field, 6' = 0. (b) B field, 9 = ^. (c): j field, 9 — 0. (d): (c): j field, 9 = ^. 



flow recirculation loops (r ~ 0.7 and z ~ ±0.5 see Fig.|21) 
we note a low level of magnetic field: it is expelled from 
the vortices. Close to the outer boundary, we mainly ob- 
serve a strong transverse dipolar field (Fig. llSlupper-left') 
correlated with two small loops of very strong current 
density j fPig. 1151 lower-left'l. These current loops seem 
constrained by the boundary, and might dissipate great 
amount of energy by Joule effect (see discussion below). 

B. Effects of the conducting layer 

As indicated in the first section, the main effect of 
adding a conducting layer is to strongly reduce the 



threshold. In Fig. ^1 we plot the critical magnetic 
Reynolds number for increasing values of the layer thick- 
ness. The reduction is important: the threshold is al- 
ready divided by 4 for w = 0.4 and the effects tends 
to saturate exponentially with a characteristic thickness 
w = 0.14 (fit in Fig. Ilti|) . as observed for an a^-model 
of the Karlsruhe dynamo by Avalos et al. ■ Adding 
the layer also modifies the spatial structure of the neutral 
mode: isodensity surface for w = 0.6 is plotted in Fig. 
nn with the corresponding sections of B and j fields in 
Fig. ^] The two twisted bananas of axial field are still 
present in the core, but the sheets of magnetic energy 
near the r = 1 boundary strongly develop. Instead of 
thin folded sheets on both sides of the equatorial plane. 
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FIG. 16: Critical magnetic Reynolds number vs layer thick- 
ness w. TM73 velocity field. Fit: R%,{w) = 38 + 
58 exp{--^) for w > 0.08. 




FIG. 17: Isodensity surface of magnetic energy (50% of the 
maximum) for the neutral mode with w — 0.6. 



the structures unfold and grow in the axial and azimuthal 
directions to occupy a wider volume and extend on both 
sides of the flow/conducting-layer boundary r = 1. This 
effect is spectacular and occurs even for low values of w. 

Small conducting layers are a challenge for numerical 
calculations: as far as the measured tangential velocity 
at the wall is not zero, adding a layer of conductor at 



rest gives rise to a strong velocity shear, which in practice 
needs at least 10 grid points to be represented. The max- 
imal grid width used is 0.005: the minimal non-zero w is 
thus w = 0.05. The exponential fit in Fig. ^| is relevant 
for w; > 0.1. We can wonder if the departure from expo- 
nential behavior is of numerical origin, or corresponds to 
a cross-over between different dynamo processes. 

The analysis of B and j profiles in Fig. 1 181 first reveals 
smoother B-lines and much more homogeneous a repar- 
tition for the current density. The azimuthal current 
loops responsible for the transverse dipolar magnetic field 
now develop in a wider space fFig. 1181 lower- left). Two 
poloidal current loops appear in this plane, closing in 
the conducting shell. These loops are responsible for the 
growth of the azimuthal magnetic field at r = 1 (Fig. IT^ 
upper-left). Changes in the transverse plane (6* =_?) 
are less spectacular. As already stated in Refs. pa. l43j . 
the positive effect of adding a layer of stationary conduc- 
tor may reside in the subtle balance between magnetic 
energy production and Ohmic dissipation. 



C. Energy balance 

In order to better characterize which processes lead to 
dynamo action in a von Karman fiow, we will now look 
at the energy balance equation: let us first separate the 
whole space into three domains. 

• Qi : < r < 1 (inner fiow domain) 

• flo : 1 < r < 1 + w (outer stationary conducting 
layer) 



> 1 + w (external insulating medium) 



In any conducting domain fl^, we write the energy 
balance equation: 



d_ 

dt 



f B^ = Rrn / (j X B).V- / f+ [ (B X E).n 



The term in the left part of the equation is the tem- 
poral variation of the magnetic energy Emag- The first 
term in the right part of the equation corresponds to 
the source term which writes as a work of the Lorentz 
force. It exists only in 17^ and is denoted W . The second 
term is the Ohmic dissipation D, and the last term is the 
Foynting vector fiux P which vanishes at infinite r. 

We have checked our computations by reproducing the 
results of Kaiser and Tilgner ^i^l on the Ponomarenko 
fiow. 

At the dynamo threshold, integration over the whole 
space gives 

Q^W-Do-D, 

In Fig. ^1 we plot the integrands of W and D at 
the threshold for dynamo action, normalized by the total 
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FIG. 18: Meridian sections of B and j fields for the neutral mode with w — 0.6. B is divided by the total magnetic energy. 
Arrows correspond to components lying in the cut plane, and color code to the component transverse to the cut plane. A unit 
arrow is set into each figure lower left corner, (a): B field, 6* = 0. (b) B field, ^ = f. (c): j field, 9 — 0. (d): (c): j field, ^ = f . 
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FIG. 19: (a): radial profile of Ohmic dissipation integrated 
other 6 and z: Jq^ J^Qgi" j^i^) dz d6 for increasing values 
of w. (b): radial profile of magnetic energy production inte- 
grated other 6 and z: Jq^ jj^^r ((j x B).V)(r) dz for 
increasing values of w. 



FIG. 20: Ratio of the integrated dissipation in the outer 
region and in the inner region vs w. Fit: -^(w) = 
0.16(l-exp(-^)). 



instantaneous magnetic energy, as a function of radius 
r for various w. For w = 0, both the production and 
dissipation mostly take place near the wall between flow 
and the insulating medium (r = 1), which could not have 
been guessed from the cuts of j and B in figure the 
w = curve in Fig. 1191 has two bumps. The first one 
at r ~ 0.1 corresponds to the twisted bananas, while the 
second is bigger and is localized near the flow bound- 
ary r — 1. A lot of current should be dissipated at the 
conductor-insulator interface due to the "frustration" of 
the transverse dipole. This can explain the huge effect of 
adding a conducting layer at this interface: the "strain 
concentration" is released when a conducting medium is 
added. So if we increase w, the remaining current con- 
centration at r = \ + w decreases very rapidly to zero, 
which explains the saturation of the effect. In the mean 
time, the curves collapse on a single smooth curve, both 
for the dissipation and the production (solid black curves 
in Fig. [Tn|l . For greater values of w, the production den- 
sity and the dissipation in the core of the flow r < 0.2 are 
smaller, whereas a peak of production and dissipation is 
still visible at the flow-conducting shell interface r = 1. 
The conducting layer does not spread but reinforces the 
localization of the dynamo process at this interface. This 
can help us to understand the process which rises the dy- 
namo in a von Karman type flow. 

Let us now look at the repartition between the dissi- 
pation integrated over the fiow Di and the dissipation 
integrated over the conducting shell Do (Fig. I20() . The 
ratio Do/Di increases monotonically with w and then 
saturates to 0.16. This ratio remains small, which con- 
firms the results of Avalos et al. 01 for a stationary dy- 
namo. We conclude that the conducting layer existence 
— allowing currents to flow — happens to be more impor- 
tant than the relative amount of Joule energy dissipated 
in this layer. 
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D. Neutral mode structure 

From the numerical results presented above in this sec- 
tion, we consider the following questions : Is it possi- 
ble to identify typical structures in the eigenmode of the 
von Karman dynamo ? If yes, do these structure play 
a role in the dynamo mechanism ? We have observed 
magnetic structures in the shape of bananas and sheets 
(see Figs. [TH and [T7|) . In the center of the flow volume, 
there is an hyperbolic stagnation point equivalent to "a- 
type" stagnation points in ABC-flows (with equal coeffi- 
cients) 13 • III the equatorial plane at the boundary the 
merging of the poloidal cells remembers "/?-type" stag- 
nation points in ABC-flows. In such flows, the magnetic 
field is organized into cigars along the a-type stagnation 
points and sheets on both sides of the /3-type stagna- 
tion points |45| |: this is very similar to the structure of 
the neutral mode we get for w = (Fig. I14|l . We also 
performed magnetic induction simulations with an im- 
posed axial field for the poloidal part of the flow alone. 
We obtain a strong axial stretching: the central stag- 
nation point could be responsible for the growth of the 
bananas/cigars, which are twisted by the axial differen- 
tial rotation after. One should nevertheless not forget 
that real instantaneous flows are highly turbulent, and 
that such peculiar stagnation points of the mean flow are 
especially sensitive to fluctuations. 

The presence of the conducting layer introduces new 
structures in the neutral mode (see Figs. 1141 1171 and [T51 
I18|l . In order to complete our view of the fields in the 
conducting layer, we plot them on the r ~ 1 cylinder for 
w = 0.6 fFig. l21|l . As for w = 0, the dipolar main part of 
the magnetic field gets radially into the flow volume at 
6 = TT and exits at 9 = (Fig. [2] up). However, looking 
around z = 0, we observe that a part of this magnetic 
flux is azimuthally diverted in the conducting shell along 
the flow boundary. This effect does not exist without 
conducting shell: the outer part of the dipole is anchored 
in the stationary conducting layer. 

Another specific feature is the anti-colinearity of the 
current density j with B at (z = 0; = 0,7r; r = 1), which 
could remind an "a" -effect. However, while the radial 
magnetic field is clearly due to a current loop (arrows 
in the center of Fig. 1211 down), jr is not linked to a B- 
loop (Fig. 1^ up), which is not obvious from Fig. [T51 
Thus, the anti-colinearity is restricted to single points 
(z = 0; 6* = 0, tt; r = 1). We have checked this, computing 
the angle between j and B: the isocontours of this angle 
are very complex and the peculiar values corresponding 
to colinearity or anti-colinearity are indeed restricted to 
single points. 

E. Dynamo threshold reduction factor 

We have shown that the threshold for dynamo action is 
divided by four when adding a conducting layer of thick- 
ness w = 0.4. This effect is very strong. Following Avalos 




FIG. 21: (a): (resp. (b)) B (resp. j) field at r = 1 for w = 0.6. 
Color code corresponds to Br (resp. jr) and arrows to Bz and 
Be (resp. and je). 

and Plunian |4^ , let us compare the threshold reduction 
factor A = 1 — JI'^^^q'^ for various kinematic dynamos. 
The threshold reduction for TM73-flow (A = 0.78) is 
much higher than for the Karlsruhe (A = 0.11) and Riga 
(A = 0.56) dynamos. Reduction rate can also be rad- 
ically different between model flows: the a^-model for 
Karlsruhe dynamo gives a low- dynamo for w = 
and benefits very few of a finite w (A — 0.11), while 
the Ponomarenko flow does not lead to dynamo action 
without conducting layer (A = 1). The reduction factors 
considered above are maximal values obtained either for 
high w in stationary dynamos or for the optimal w in 
oscillatory dynamos |42. |43| . 

In order to understand why A is so high for our TM73- 
flow, we propose to compare our experimental flow with 
an optimal analytical model-flow proposed by Marie, 
Normand and Daviaud ^39] in the same geometry. The 
Galerkin method used by these authors does not allow to 
study the effect of a conducting layer. We thus perform 
kinematic dynamo simulations with our usual approach, 
and then study the effects of adding a conducting layer on 
the following velocity field for e — 0.7259 corresponding 
to r = 0.8 [^IH: 

Vr = -^r(l - r)^(l 2r) cos(7r2;) 

vg — 4er(l — r) sin(7rz/2) 

= (1 - r)(l + r - 5r^)sin(7rz) (2) 

The kinematic dynamo threshold is found at i?^ = 58 
for w = 0, in good agreement with the galerkin analysis. 
With a w — 1 conducting layer, we get a low A — 0.26 
reduction rate, i.e. = 43, close to the TM73 thresh- 
old ioT w — 1: = 37. The threshold reduction is 
also found to show an exponential behavior with w, of 
characteristic thickness 0.20, as in Fig. 1161 

Let us describe the model flow features represented 
in Fig. [7| (bottom). The velocity is very smooth at the 
cylindrical boundary: the toroidal velocity is maximum 
at r = 0.5 and slowly decreases to zero at r = 1. The 
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poloidal recirculation loops are centered at = 0.56 
and the axial velocity also decreases slowly to zero at 
the cylindrical boundary. Thus, the mass conservation 
requires the axial velocity to be much higher in the cen- 
tral disk (0 < r < Vp) than outside. These constraints 
make analytical models somewhat different from exper- 
imental mean flows (Fig. 0)- In particular, high kinetic 
Reynolds numbers forbid smooth velocity decrease near 
boundaries. This explains why experimental flows do not 
lead to low thresholds unless a conducting layer is added. 

We now consider the effect of a conducting shell on the 
model flow eigenmode structure. First note that with- 
out conducting shell, the model neutral mode structure 
is already very similar to that of TM73 with conduct- 
ing shell: the transverse dipole is not confined into thin 
sheets but develops into wider regions connected to ba- 
nanas of axial field in the center. Adding the conducting 
layer mainly lets the neutral mode structure unchanged 
and thus quantitatively reduces its impact compared to 
the experimental case. 

Finally, from the very numerous simulations of experi- 
mental and model von Karman flows performed, we con- 
clude that the adjunction of a static conducting layer to 
experimental flows makes the eigenmode geometry closer 
to optimal model eigcnmodes, and critical i?^ get closer 
to moderate values (typically 50). It may thus be conjec- 
tured that the puzzling sensitivity of dynamo threshold 
to flow geometry is lowered when a static layer is present. 
We conclude this feature renders the dynamo more robust 
to flow topology details. This could also act favorably in 
the nonlinear regime. 



V. CONJECTURES ABOUT DYNAMO 
MECHANISMS 



present azimuthal velocity shear in axial and radial di- 
rections (see Fig.|21l. So, we will consider below the role 
of both axial and radial w-effect. 

We will discuss these mechanisms and then suggest 
that, for a flow surrounded by a static conducting layer, 
the dynamo mechanism is based on the presence of a 
strong velocity shear (at the boundary layer r = 1) which 
lies in this case in the bulk of the overall electrically con- 
ducting domain. 

A. Axial a;-efFect 

Induction simulations performed with the toroidal part 
of the velocity show an axial w-effect which converts an 
imposed axial field into toroidal field through dve/dz. 
Such a i?Tp,-linear effect has been evidenced in VKSl ex- 
periment [2y|. This effect concentrates around the equa- 
torial shear layer {z = 0) as visible in Fig. |21 Thus, we 
can think that the axial w-effect is involved in the dy- 
namo process: for dynamo action to take place, there is 
a need for another process to convert toroidal magnetic 
field into poloidal field. 

B. Q-efFect, helicity effect 

i?„i-iion-linear conversion from transverse to axial 
magnetic field has also been reported in VKSl experi- 
ment |2]J. This effect is not the usual scale-separation 
Qf-effect |2g| and has been interpreted as an effect of the 
global helicity as reported by Parker (in the follow- 
ing, it will be denoted "a" -effect). We believe it to take 
place in the high kinetic helicity regions of the flow (see 
Fig.lHl). 



In this paragraph, we intend to relate the results of 
the optimization process to some more elementary mech- 
anisms. As emphasized in the Introduction, there is no 
sufficient condition for dynamo action and although nu- 
merical examples of dynamo flows are numerous, little is 
known about the effective parameters leading to an effi- 
cient energy conversion process. For example, the clas- 
sical a and axial lo mechanisms have been proposed to 
be the main ingredients of the von Karman dynamo |l9j | . 
Our starting point is the observation that dynamo ac- 
tion results from a constructive coupling between mag- 
netic fields components due to velocity gradients, which, 
in the present axisymmetric case, reduce to derivatives 
with respect to r (radial gradients) and to z (axial gra- 
dients). The gradients of azimuthal velocity generate a 
toroidal field from a poloidal one ((jj-effect JJ), while re- 
generation of the poloidal field is generally described as 
resulting from an helicity eflFect (denoted a-effect if scale 
separation is present [2g). How do these general con- 
siderations apply to the present flow ? As in the Sun, 
which shows both a polar-equatorial differential rotation 
and a tachocline transition, our experimental flow fields 



C. Is an "a"[i; mechanism relevant ? 

Bourgoin et al. |46j | performed a study of induction 
mechanisms in von Karman-type flows, using a quasi- 
static iterative approach. They show that "q;"w dynamo 
action, seen as a three-step loop-back inductive mecha- 
nism, is possible, but very difficult to obtain, fields being 
widely expelled by the vortices. The authors highlight 
the fact that the coupling between the axial w-effect and 
the "a" -effect is very inefficient for our velocity fields, be- 
cause of the spatial separation of these two induction ef- 
fects. Our observations of the velocity and helicity fields 
confirm this separation. 

The authors also discovered an induction effect — the 
BC-effect — related to the magnetic diffusivity disconti- 
nuity at the insulating boundary that could be invoked 
in the dynamo mechanism. This BC-effect. illustrated on 
our TM73- velocity field (Fig. 14 in Ref. [ig), is enhanced 
in the case of strong velocity and vorticity gradients at 
the boundaries, characteristic of high Reynolds number 
flows. So, we are convinced that for experimental flow 
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fields at w = 0, the BC-effect helps the dynamo. This is 
coherent with our observations of high tangential current 
density near the boundaries and high magnetic energy 
production at r = 1 even for w = (Fie. Such a 

current sheet formation and BC-effect was reported by 
BuUard and Gubbins 

When a large layer of sodium at rest is added, the 
BC-effect vanishes because the conductivity discontinu- 
ity occurs at r = 1-l-w while the currents still concentrate 
at the flow boundary r = 1. However, with a conduct- 
ing layer, we have presented many features favoring the 
dynamo: in the next paragraph, we propose a possible 
origin for this conducting-layer effect. 

D. Radial w-effect, boundary layers and static shell 

With a layer of steady conducting material surround- 
ing the flow, we note the occurrence of two major phe- 
nomena: 

• the possibility for currents to flow freely in this shell 

• the presence of a very strong velocity shear local- 
ized at the boundary layer which now lies in the 
bulk of the electrically conducting domain. 

Let us again consider the shape of the velocity shear. 
Any realistic (with real hydrodynamical boundary condi- 
tions) von Karman flow obviously presents negative gra- 
dients of azimuthal velocity dve jdr between the region 
of maximal velocity and the flow boundary. This region 
can be divided into two parts: a smooth decrease in the 
bulk {R r < 1) and a sharp gradient in the boundary 
layer at r = 1 (Fig. [7|). 

These gradients are responsible for a radial w-effect, 
producing Bg with Br, in both insulating and conduct- 
ing cases. However, without conducting layer, only the 
smooth part of the gradient which lies in the bulk will be 
efficient for dynamo action. Indeed, owing to the huge 
value of the kinetic Reynolds number and the very small 
value of the magnetic Prandtl number, the sharp bound- 
ary layer gradient is confined in a tiny domain, much 
smaller than the magnetic variation scale. No significant 
electrical currents can flow in it and we did not resolve 
this boundary layer with the numerical code: it is totally 
neglected by our approach. 

The role of both types of gradients is illustrated by 
the observation (Fig. [7| left) of impellers of large radius 
(i? — 0.925). For such impellers there is almost no de- 
parture from solid body rotation profiles in the flow re- 
gion and these impellers lead to dynamo action only with 
conducting shell 22], i.e., due to the sharp gradient. On 
the other hand, our R — 0.75 selected impellers present a 
stronger bulk-gradient and achieve dynamo in both cases. 

Actually, the way we numerically modelized the von 
Karman flow surrounded by a static conducting layer 
— considering an equivalent fluid system in which the 



boundary layer appears as a simple velocity jump in its 
bulk — is coherent with the problem to solve. The veloc- 
ity jump, just as any strong shear, is a possible efficient 
source for the radial tj-effect. 



E. A shear and shell dynamo ? 

We pointed out above that the regions of maximal he- 
licity (the "a" -effect sources, see Fig. I12|) are close to 
those of radial shear where radial w-effect source term 
takes place. Dynamo mechanism could thus be the result 
of this interaction. So, in the absence of a static shell, 
one can suppose that the dynamo arises from the cou- 
pling of "a" -effect, w-effect and the BC-effect [i^. With 
a static conducting layer, as explained above, the radial 
w-effect is especially strong: the radial dipole, anchored 
in the conducting layer and azimuthally stretched by the 
toroidal flow (see Fig. l21|) is a strong source of azimuthal 
fleld. This effect coupled with the "a" -effect could be at 
the origin of the dynamo. 

For small conducting layer thickness w, one could ex- 
pect a cross-over between these two mechanisms. In fact, 
it appears that the decrease of R'^ (Fig.ll6|l with the con- 
ducting shell thickness w is very fast between w — and 
w = 0.08 and is well fitted for greater w by an exponen- 
tial, as in Ref. We can also note that for typical 
Rm — 50, the dimensionless magnetic diffusion length 

— 1/2 

Rm is equal to 0.14. This value corresponds to the 
characteristic length of the decrease (Fig. and is 
also close to the cross-over thickness and characteristic 
lengths of the Ohmic dissipation profiles (Figs. ^1 top 
andHni)- 

We propose to call the mechanism described above a 
"shear and shell" dynamo. This interpretation could also 
apply to the Ponomarenko screw-flow dynamo which also 
merely relies on the presence of an external conducting 
medium. 



VI. CONCLUSION 

We have selected a configuration for the mean flow fea- 
sible in the VKS2 liquid sodium experiment. This mean 
flow leads to kinematic dynamo action for a critical mag- 
netic Reynolds number below the maximum achievable 
Rm- We have performed a study of the relations be- 
tween kinematic dynamo action, mean flow features and 
boundary conditions in a von Karman-type flow. 

The first concluding remark is that while the dynamo 
without static conducting shell strongly depends on the 
bulk flow details, adding a stationary layer makes the 
dynamo threshold more robust. The study of induction 
mechanisms in 3D cellular von Karman type flows per- 
formed by Bourgoin et al. ji^ suggests that this sensi- 
tivity comes from the spatial separation of the different 
induction mechanisms involved in the dynamo process: 
the loop-back between these effects cannot overcome the 
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expulsion of magnetic flux by eddies if the coupling is 
not sufficient. Secondly, the role of the static layer is 
generally presented as a possibility for currents to flow 
more freely. But, instead of a spreading of the currents, 
the localization at the boundary of both magnetic energy 
production and dissipation (Fig. I19|l appears strongly re- 
inforced. Actually, strong shears in the bulk of the elec- 
trically conducting domain imposed by material bound- 
aries are the dominating sources of dynamo action. They 
result in a better coupling between the inductive mech- 
anisms. We also notice that there seems to be a general 
value for the minimal dynamo threshold (typically 50) 
in our class of flows, for both best analytical flows and 
experimental flows with static conducting layer. 

Although the lowering of the critical magnetic 
Reynolds number due to an external static envelope 
seems to confirm previous analogous results 0, ' 
it must not be considered as the standard and general 
answer. In fact, in collaboration with Frank Stefani 
and Mingtian Xu from the Dresden MHD group, we are 
presently examining how such layers, when situated at 



both flat ends, i.e., besides the propellers, may lead to 
some increase of the critical magnetic Reynolds number. 
This option should clearly be avoided to optimize fluid 
dynamos similar to VKS2 configuration. However, a spe- 
cific study of this later effect may help to better under- 
stand how dynamo action, which is a global result, relies 
also on the mutual effects of separated spatial domains 
with different induction properties. 
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